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ABSTRACT 

The correlation between radio and X-ray luminosity in the hard state of black hole X- 
ray binaries is important for unveiling the relation between the accretion flow and the 
jets. In this paper, we have modeled the quasi-simultaneous multi-band observations 
of a recently discovered transient X-ray source, IGR J17177-3656. It is found that the 
source is probably an outlier following the steep radio/X-ray correlation rather than an 
outlier in the transition region as suggested by Paizis et al. (2011). It is also found that 
the multi-band spectrum can be successfully modeled by the luminous hot accretion 
flow (LHAF) but less likely by the advection dominated accretion flow (ADAF). Our 
results support the point that LHAF can explain the steep radio/X-ray correlation. 

Key words: X-rays: binaries - X-rays: individual: IGR J17177-3656 - stars: winds, 
outflows - black hole physics - accretion, accretion disks 



Fender I l2006l : [Pone et al\ 12003). During the hard 
there exists an interesting nonlinear correlation be- 



1 INTRODUCTION 

The hard state of the black hole (BH) X-ray binaries is 
characterized by a power-law X-ray spectrum with photon 
inde x F = 1.5 — 2.1, a nd relatively stro ng radio emission 
(e.g.lZdziarski fc Gierlhi ski 2004; Rcmilla rd fc McClintodTI 
l2006l : iFender I l2006l : [Pone et al\ 120071 ). 
states 

tween the radio and X-ray luminosity, i.e. Lra d jo oc -t/y with 
b ^ 0.5 - 0.7 (e.g. iHannikainen et al\ Il998l : ICorbel et all 
I2OOOI . I2OO3I : iGallo et aZ.ll2003l ). Although this relation is es- 
tablished from a few BH binaries like GX 3 39-4 and V404 
Cyg, i t even holds for ac t ive galactic nuclei IJMerloni et al\ 
I2OO3I : iFalcke et al} l2004l : iGiiltekin et~al\ |2009| ). Moreover. 
a similar correlatio n between o ptical /infrare d and X-ray 
flux has been foun d (jHoman et a l. 2005; Russe ll et fflZ.ll2006l : 
ICoriat et aZ.ll2009l '). It seems that this correlation is universal 
and can be taken as "standard" . 

However, radio observations of some new BH tran- 
sients in more recent outbursts in the following years reveal 
that their behavior devi ates significantly from the aforemen- 
tioned correlation (e.g . I Corbel et al. 20041 : [Rodriguez et al\ 
l2007l : IXue fc Cui ll2007l : [Soleri et oZ.l2010l ). Bas ed on the first 
precis e measurement of an outlier H1743-322, ICoriat et al\ 
(|201ll . Cll) discovered a new tight correlation in its bright 
hard state, with the slope being 1.38 ±0.03. At low luminos- 
ity H1743-322 rejoins the standard correlation with b ~ 0.6. 
At the intermediate luminosity, there exists a transition re- 
gion. The reasons for the existence of the two correlation 
and their transition remain to be clarified. 
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To date three classes of models have been pro- 
posed to explain the hard states: (i) the hot accretion 
flow, which includes ADA F (e.g. 
Abramowicz et al. 1995; Chen et al. 



Yuan | |2007|; iNaravan fc McCUntock 



Yuan I I2OO 
200411 ■ (ii) 



J, I2OO3I ; lYuan et al 
the dis k coro n a (e.g 



Naravan fc Yi I 



1995; Na ravan I 



1994; 



20051 : 



20081') and LHAF (e.g 



20qi: lYuan fc ZdziarskTI 
Liang fc Pricel 



1977 



2OOOI : iLiu e^ al\ 



Haardt fc Maraschi 1991; Zhang et al\ 
2002, 2007, 2011; Rcis et al. 2010) , and ^. 

2001, 20 0A iReig eTai\m m-. Gi,^ 

2009i r |Pe'er fc Casella .,20091 ). All these 



(iMarkofi' et al.\ 
I2OO4I : 



Maitra et 



(iii) the jets 
et al\ 



models can explain the "standard" correlation between ra- 
dio and X-ray emission. In the models of hot accretion flow 
and corona. X-rays come from the accretion flow while radio 
emission from the jet. Since the jet originates from the ac- 
cretion fl ow, it is natural to expect the radio/X-ray correla- 
tion (e.g. [Merloni et fflZ.l2003l : lYuan fc Cui l2005l : lYuan et al] 
'2009'). In the jet model, both radio and X-ray emission come 
from the same elect rons in the jet, and thus strong correla - 
tion is expected (e.g. lMarkoff ei aZ.l200l : lFalcke et d.l2004l) 
As for the steep correlation with index of ~ 1.4. ICoriat et al\ 
(|2011l ) proposed that the efficient accretion flow like 
LHAF or disk corona can explain it. To understand 
the physics behind the correlation, more sample of sources 
are needed. 

Very recently, I Paizis et al\ (|201ll ) (Pll) reported the 
observation of a newly discovered X-ray transient, IGR 
J17177-3656. They analyzed the quasi-simultaneous data 
by Chandra, INTEGRAL and ATCA, and argued that the 
source is probably a low-mass BH X-ray binary in the hard 
state, although the neutron star nature cannot be excluded 
yet. By assuming the distance to be 8 kpc, they found that 
the source follows the behavior of the BHC H1743-322 at 
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intermediate X-ray luminosity. In this paper, by using an 
accretion-jet model, we obtain the source distance by model- 
ing the observed multi-band spectrum, and re-investigate its 
position in the radio/X-ray correlation diagram. It is found 
that IGR J17177-3656 probably follows the steep radio/X- 
ray correlation of H1743-322 and not the transition one. 



2 MODELING THE SPECTRUM 

2.1 Observational data 

After the discovery by INTEGRAL, IGR J17177-3656 be- 
came a focus of radio, NIR, and X-ray telescopes (Pll 
and references therein). Chandra conducted an observa- 
tion of 20 ks at March 22nd, 2011. The 2-8keV observa- 
tion shows a power-law spectrum, F — 1.36lo'i5i with un- 
absorbed 3-9 keV flux being 7.5 x 10"^"^ erg cm~^ s~^. A 
simultaneous INTEGRAL observation in 20-200 keV shows 



S -1-0. 94 



a power law with V — 1.8_i'4 and flux being QA8_{ 
10~^° erg cm~'^ s~^. This result is similar to the observa- 
tion on March 15th, 2011. The radio observation on the 
same day is conducted by Australia Telescope Compact Ar- 
ray (ACTA) at two frequency bands, which are centered at 
5.5 GHz and 9 GHz. The corresponding flux density are 
0.24 ± 0.06 mjy and 0.20 ± 0.06 mjy, respectively. Figure 1 
shows the results of our modeling and the data points of 
Pll. 



2.2 The accretion-jet model 

The accreti on-jet model is br iefly described here. The reader 
can refer to lYuan et al\ ((2005,) for additional details. There 
are four components in the model, the inner hot accretion 
flow (ADAF/LHAF), the outer standard accretion disk, the 
outflow and the jet. In the ADAF regime, the accretion 
rate M < a^Msdci, where MEdd ~ WLEdd/c^ and a is 
the viscous parameter, and the viscously dissipated energy 
can be advected inwards with the accreting gas. While in 
the LHAF regime, a^MEdd < M < 5a^MEdd, and the 
radiative cooling is a little more efficient than the viscous 
heating. Although the temperature of LHAF is lower than 
that of ADAF, it is also high because during the falling 
the compression work plus viscous heating together balance 
the radiative cooling. The jet is shown in th e flat /inverted 
radio spectra and the di r ect images (e . g. Fender I l200ll . 
l2006l : iDhawan et~ai\ l2000l: IStirUng et ail 1200 ll ). The exis- 
tence of outflow in ADAF/L HAF is from theoretical work or 
numerical simulations (e.g. Ilgumenshchev fc Abramowicz I 
1999"; 'Naravan & Yi ' 'l994'; ^Blandfo rd fc Begelman 1 1 19991 : 
Narayan et al. 2000; Yuan & Bu 201^). The outer standard 
accretion disk is needed w hen modeling the UV spectra of 
some well-known BHs (e.g . iNaravan et ffl/.lll996l ; lEsin et ali 
I2OOII : iFrontera ei"aZll200ll ) . 

The radiation processes in the hot accretion flow include 
synchrotron, bremsstrahlung, and inverse Compton scatter- 
ing of soft photons. The emission from the outer cool disk 
is modeled as a multi-color blackbody spectrum. And in the 
jet, synchrotron emission and self-absorption of the relativis- 
tic electrons are considered. 

The accretion disk is described by a set of parameters. 
The inner hot accretion flow truncates the outer standard 



accretion disk at the transition radius rtr- Since the physics 
of the transition is highly uncertain, the range of the value of 
rtr i s large, from ^ l Org to ~ 10'* rg, where rs ~ 2GM/c? 
(e.g. lEsin et aZ.1 119971 ). Due to the lack of UV spectra, we 
set it a normal value, rtr = lOOrg. The influence of out- 
flow on the accretion is mimicked with the parameter s by 
M = Mo(r/rtrY , where Mq is the accretion rate at rtr- The 
value of s describes the strength of the outflow and varies 
in the range < s < 0.5. For the case s = 0, there is no 
outflow, while for the case s — 0.5, 90% of the accreted gas 
will flow away and only 10% flnally fall into the BH given 
rtr = lOO rg, which roughly agre es with numerical simula- 
tions (e.g. lStone fc Pringle IbOOll ). A recent numerical simu- 
lation shows that the f raction of outflow is almost the same 
for ADAF and LHAF (JYuan fc Bull2010l '). So in this paper, 
the value of s is assumed to be the same for different accre- 
tion rates, i nstead of a f unctio n related to accretion rate as 
assumed in lYuan et al\ (120051 ). To determine the dynamics 
of the hot accretion flow, we need three more parameters, 
which are associated with the MHD turbulence. These pa- 
rameters are the viscous parameter a, the fraction of gas 
pressure to total pressure /3, the fraction of the viscous dis- 
sipation directly heating electrons S. The value of /3 is fixed 
at a typical value of 0.9, the v alue of a changes in the rang e 
0.01 - 0.3, and S ~ 0.01 - 0.5 (|Q.uataert fc Gruzinov Ill999l ). 

The jet is modeled phenomenologically with quite a 
few parameters. We assume a fraction of the accretion flow, 
rhjet, is transferred into the vertical direction, and the jet 
has a conical geometry with a small half-opening angle ?/). 
For simplicity we assume each ejection moves at a constant 
velocity along the jet. Different ejections may have differ- 
ent velocities, however, their averaged bulk Lorentz factor, 
Fj, is assumed to be constant. When different ejections col- 
lide, internal shocks occur, which accelerate a fraction of 
the electrons into power-law energy distribution with index 
p. The energy density of accelerated electrons and amplified 
magnetic field are specified with ee and es. lYuan fc Cui I 
(2005) have shown that the contribution of the jets to the 
X-ray emission in BH sources is negligible when the lumi- 
nosity is > 10~^ LEdd- Therefore in the transient source IGR 
J17177-3656, the jet contributes only to the radio emission. 
It should be noted that the main results of the paper are 
from the modeling of the X-ray observation, and we model 
the radio spectra just to show the ability of the model to 
explain the multi-band spectra. So we fix most of the jet 
parameters. Considering that there are only two radio data 
points, we leave rhjet and p free and the other parameters 
fixed at the value that has been used for other BH binaries, 
i.e. -0 = 0.1, F, = 1-2, e e = 0.06 and es = 0.02 (|Yuan et al\ 
I2OO5I : IZhang et aLll2010l ) . 

The BH mass is set it to be of normal value A'Ibh = 
IOM0 . The binar y inclination angle 8 is estimated to be high 
60° < e < 75° (jPaizis et a^.ll2011^ . Detailed calculations 
show the influence of 9 is not significant, so we set it to be 
a middle value d = 70°. 

There are 7 free parameters and quantities in the model, 
Af, s, a, S, TTijet, p, and the distance d. The modeling steps 
are as follows. First, we try with different value of M, s, a, 
and S to model the slopes of the X-ray spectra, and obtain 
the possible ranges of these parameters. Second, the range 
of d is obtained by comparing the observed X-ray flux and 
theoretically calculated disk luminosity. At last, nijet and 
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Figure 1. Spectral modeling results for the quasi-simultaneous 
observation of IGR J17177-3656 by ATCA (5GHz, 9GHz), 
C/iara(ira(3-9 keV) and INTEGRAL (20-200 keV). The thick solid 
line is the total spectrum and corresponds to Case I with s = 0.2, 
M = 0.4M£;rfrf, and d = 25 kpc. For Case I, the short dash-long 
dashed, long dashed, and dot-dashed lines correspond to radiation 
from the jet, the outer thin disk, and the inner hot accretion flow, 
respectively. The thick dashed line is the total spectrum and cor- 
responds to Case II with s = 0.3, M = 0.5M^^^, and d = 21 kpc. 
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Figure 2. The variation of the advection factor. The solid and 
dashed lines correspond to Case I and II, respectively. 



p are adjusted to check whether the radio spectrum can be 
modeled simultaneously. 



2.3 Modeling results 

It is useful to investigate the effects of parameters on the 
output spectra in order to find the possible ranges of the pa- 
rameters. Note that the X-ray spectral index decreases with 
the increasing y Compton parameter, which is the product 
of electron temperature and optical depth. Given other pa- 
rameters fixed, the surface density of the accretion fiow of 
higher M is higher, therefore the optical depth is higher 
and the spectral index is smaller. The value of a determines 
the maximum accretion rates of the ADAF and LHAF, and 
therefore the maximum optical depth and hardest spectra. 
For larger 5, more viscously dissipated energy is transferred 
to thermal electrons. Consequently, the electron tempera- 



Figure 3. Correlation between radio and X-ray luminosity for a 
sample of black holes and neutron stars. The asterisk and dotted 
line mark the possible location of IGR J17177-3656. The other 
data are adapted from Cll. The solid and dashed lines are 
the fit with correlation index of 0.6 and 1.4. 



ture is a little higher and the X-ray spectrum is harder. 
When s increases, the final X-ray spectrum is softer even if 
the electron temperature increases due to less cooling. This 
is because, as shown by detailed calculations, the decrease 
of the surface density and optical depth are more significant 
than the increase of the electron temperature. It should be 
mentioned that if the X-ray luminosity rather than M is 
given, the X-ray spectral hardness increases with s. This is 
because more photons come from t he region of larger radius, 
where the optical depth is higher (lYuan et fflZ.ll2007l ). 

Our detailed calculations show that in order to model 
the observed spectrum, the parameters are limited in the 
following ranges, M > O.SMEdd, ce > 0.2, 5 > 0.1, and 
0-2 < s < 0.3. When M < 0.3MEdd, a sj 0.2 or 5 < 0.1, the 
output photon index are greater than 1.54. When s < 0.2, 
the obtained luminosity is as large as 8.1 x lO^'^ erg s~^, and 
the obtained distance is greater than 30 kpc, which seems 
unreasonable since the diameter of the Galaxy is about 
30kpc. When s > 0.3, the highest photon index obtained 
with M = bo?MEdd is 1.65 , which is not hard enough to 
explain the observed X-ray spectra. 

The modeling is not unique. As an example two cases 
of the modeled spectra are shown in Figure 1. In Case I, 
the set of parameters are s = 0.2, a = 0.3, 5 = 0.5, M = 
QAMsdd, nijet = 0.3%, p = 2.25, and d = 25 kpc. While 
in Case II, s = 0.3, M = O.bMEdd, and d = 21 kpc, with 
others being the same as Case I. The thick solid and dashed 
lines correspond to Case I and II, respectively. From this 
figure it can be seen that the radio and X-ray emission are 
dominated by the jet and the hot accretion flow, respectively. 
The infrared and optical spectra are the sum of the emissions 
from the jet, hot accretion fiow and the truncated thin disk. 
It can also be seen that the model parameters can be further 
constrained by optical/UV observation. 

The accretion rate is in the range o?MEdd < M < 
5a^MEdd, therefore the accretion fiow is in the regime of 
LHAF. This regime is characterized by the negative advec- 
tion factor, which is defined as fadv = Qadv/q^ ~ (<?^ — 
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qie)/q^ , with (/"*", qadv, Qie being the viscous heating rate, 
advection cooling rate, and the coohng rate of ions by ion- 
electron coupling, respectively. It is easy to find in Figure 2 
that fadv < except in the innermost region r < 6rs- The 
key reason why LHAF can model the spectra but ADAF 
cannot is because of the observed hard X-ray spectra with 
photon index being ~ l.SGlgjg, which demands optically 
thick hot gas that is difficult for ADAF to form. 

The distance is 21-25 kpc, rather large compared to 
most of the observed BH X-ray binaries. However, this is 
still smaller than that of the farthest dynamically confirmed 
BH GS 1354-64, which i s estimated to lie at 25-60 kpc 
iRevnolds fc Miller llJoTH ). Therefore it is still well accept- 
able. Considering the diameter of the Galaxy is 20-30 kpc, 
we think that the source may be located on the side of the 
Galaxy opposite to the sun. Moreover, if the mass of the BH 
is adjustable, the distance could be as close as ~13-15 kpc 
for Mbh = 3Mq. 



IWrobel et fflZ.l l2008': 'Pel legrini ei"aLll2007l : lYuan et aZ.llioogI ) 
and was comp atible with theoretical works based on obser- 
vatio n al data JWu ei a/.ll2007l: Fszota e^ fflZ.ll2008l : IYu eraU 
I2OIII : iRodrigo et al\ |2012| ). iGallo et al\ l|2006h . however, 
found that the radio and X-ray data of the quiescent state 
of A 0620-00 are consistent with an extrapolation of the V 
404 Cyg and GX 339-4 correlations. Sinc e the exact track of 
A 0620-00 in out burst is still unknown jCorbel et fflZ.ll2008l : 
lYuan et ali 12003 ) , it is doubtful to say this r esult conflicts 
with the prediction of IVuan fc Gui I (|2005l '). iGorbel et~ai\ 
(2008) reported the radio and X-ray observations to the qui- 
escent state of V404 Cyg and found that the correlation of 
this source remain the original b ~ 0.6 one. This result may 
also not contradict with the 1.23 correlation, considering the 
quiescent luminosity of t his source is very c lose to the critical 
luminosity of IO^'^LecM (|Yuan et al\\200^ ). 



3 DISCUSSION 

3.1 Location of IGR J17177-3656 in radio/X-ray 
correlation diagram 

From the spectral modeling results, the distance of the 
source should be 13-30 kpc by considering all possible BH 
mass and the radius of the Galaxy. The possible position of 
IGR J17177-3656 is shown in Figure 3. It can be seen that 
except for low mass BH the possible location of the source 
is close to the branch with 6 ~ 1.4. This indicates that IGR 
J17177-3656 is a new source like H1743-322, which follows 
the steep radio/X-ray luminosity correlation. Moreover, our 
results show that LHAF can not only explain the relation 
between the radio and X-ray luminosity, but also model the 
multi-band spectrum. This supports the result of Cll that 
LHAF is a candidate to explain the correlation index with 
6~ 1.4. 

The observation of H1743-322 shows that the correla- 
tion index is ^ 0.6 at low luminosity and ~ 1.4 at high 
luminosity. The change of the correlation index can be nat- 
urally explained by virtue of the accretion-jet model. When 
A/ < a^Msdd, the hot accretion flow is in the regime 



of ADAF, the lum inosity is low and the index is 



0.6 



I Yuan fc Gui Il2005h . When M increases within a few times, 
the hot accretion flow is in the regime of LHAF, the lumi- 
nosity is higher and the index is ~ 1.4. 

LHAF is more luminous than ADAF, therefore it is used 
to explain luminous hard state with Lx > O.lLEdd- A lumi- 
nous hard state occurs during the rising part of an outburst, 
corresponding to the upper-right corner in the hardness- 
intensity diagrams (HIDs). Considering the above results it 
is natural to expect prevalent steep radio/X-ray correlation 
in the luminous hard state. Some observations seem to sup- 
port this result (Yuan, private communication). 

It should be noted that in addition to the change 
of correlation index at high luminosity, the accretion-jet 
model predicts another change of th e correlation index at 
low luminosity. lYuan fc Cui I ((20051) predicted that when 
L < 10~^LEdd, the X-ray emission should be dominated by 
the jet rather than the accretion flow. Consequently the cor- 
relation index changes from b~ 0.6 to ~ 1.23. This predic- 
tion was well confirmed by later direct observations of AGNs 



3.2 Limitations of the model 

It should be noted that the results of this paper depend on 
the observed slope of the X-ray spectrum. In the above cal- 
culation we take the upper limit of the X-ray photon index 
to be 1.52, which is of 68% confidence. If the upper limit is 
extended to 1.6, which is of 90-99% confidence as shown in 
Figure 4 of Pll, the above results are less robust. In this 
case the ADAF with weak outflow, s < 0.1, can marginally 
model the X-ray spectrum. Moreover, the minimal possi- 
ble distance of IGR J17177-3656 reduces to 18 kpc, which 
means the source locates around the point where the corre- 
lation begins to deviate from the 6 ~ 1.4 branch. 

The physics of the transition from the outer standard 
accretion disk to the inner hot accretion flow remains un- 
clear. Although different scenarios may explain the transi- 
tion, such as evaporation, turbulent dif fusive heat t ransport 
and secular instability in the cold disk (jYuan Il200ll . and ref- 
erences therein), the position and gas temperature cannot 
be estimated self-consistently. In this paper, the transition 
radius is simply assumed to be lOOrg and the ion tempera- 
ture to be 0.5-0.8 times of the Virial temperature. 

If the variation of IGR J17177-3656 is similar to H1743- 
322, there will be a transition region between the radia- 
tively inefficient and efficient branches. This means that a 
radio-quenching or X-ray-increasing process occurs during 
the transition from ADAF to LHAF. However, a physical 
explanation for such processes has not been well understood. 

In this paper, the outflow only affects the density of the 
accretion flow. However, its contribution to the loss of angu- 
lar mo mentum and energy can produce difference by a factor 
of few l|Xie fc Yuan 1120081 ). Moreover, with the existence of 
ordered magnetic field, th e loss of angula r momentum can 
be even more significant (|Bu et al\\2009^ . These effects of 
outfiow on the critical accretion rates and output spectrum 
will be studied in our future work. 

Our modeling is based on the argument of Pll that 
IGR J17177-3656 is regarded as a BH X-ray transient in 
hard state. However, the neutron star nature of the source 
cannot be excluded yet. More observations on the source are 
needed to test these results. 
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4 CONCLUSION 

In this paper we model the quasi-simultaneous multi-band 
spectrum of a newly discovered hard X-ray transient, IGR 
J17177-3656, using an accr etion-jet model. Th e model is al- 
most the same as that in lYuan et al\ (120071 ). but with a 
constant outflo w parameter s cons idering the recent numer- 
ical simulation (jYuan fc Bu Il20ld ). The required accretion 
rate is probably in the regime of LHAF, and the distance is 
about 21-25 kpc for Mbh ~ IOM0. All possible locations of 
IGR J17177-3656 in the radio/X-ray luminosity diagram are 
shown in Figure 3. This source is close to the steep correla, - 
tion branch with b ~ 1.4. The point of ICoriat et al\ ()201ll ) 
is supported, that LHAF can explain the steep radio/X-ray 
correlation. 
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